Improved calculation of relativistic shift and isotope shift in Mg I 
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We present an ab initio method of calculation of isotope shift and relativistic shift in atoms 
with a few valence electrons. It is based on an energy calculation involving combination of the 
configuration interaction method and many-body perturbation theory. This work is motivated by 
analyses of quasar absorption spectra that suggest that the fine structure constant a was smaller at 
an early epoch. Relativistic shifts are needed to measure this variation of a, while isotope shifts are 
needed to resolve systematic effects in this study. The isotope shifts can also be used to measure 
isotopic abundances in gas clouds in the early universe, which are needed to study nuclear reactions 
in stars and supernovae and test models of chemical evolution. This paper shows that isotope shift 
in magnesium can be calculated to very high precision using our new method. 
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I. INTRODUCTION 

The motivation for this work comes from recent studies 
of quasar absorption spectra designed to probe a in the 
distant past. Atomic transition frequencies depend on 
a, and by comparing frequencies on Earth with those in 
quasar absorption spectra, one can deduce whether or 
not a was different in the early universe. While some 
studies have revealed a significant deviation from zero 
(Refs. 0, 0, 0, 011other groups using a different telescope 
have not (Refs. 

One of the possible major sources of systematic effects 
in these studies is that the isotopic abundance ratios in 
gas clouds in the early universe could be very different to 
those on Earth. A "conspiracy" of several isotopic abun- 
dances may provide an alternative explanation for the 
observed variation in spectra 0, H, 0] ■ In order to test 
this possibility it is necessary to have accurate values 
for the isotope shift (IS) for the relevant atomic transi- 
tions TS]. Experimental data are available for very few 
of them; therefore, accurate calculations are needed in 
order to make the most comprehensive analysis possible. 

The need for accurate isotope shifts is further moti- 
vated by a wish to study isotopic evolution in the uni- 
verse. The isotopic abundances of gas clouds may be 
measured independently of a variation in a |10| . This is 
important for testing models of nuclear reactions in stars 
and supernovae and of the chemical evolution of the uni- 
verse. 

Previously we have calculated isotope shift in atoms 
with one valence electron using many-body perturbation 
theory (MBPT) [l^, and for neutral magnesium using 
the configuration interaction (CI) method Both of 



these papers used the finite-field scaling method for the 
perturbation. In a sense our method combines these two 
methods by including core-valence correlations into our 
CI method using MBPT. Magnesium is one of the sim- 
plest and well studied two-valence-electron atoms. Be- 
cause of that it is often used as a test ground for different 
methods of atomic calculations. In this paper we show 
that we can calculate the isotope shift of some magnesium 
transitions for which experimental values are available. 

II. METHOD 

The isotope shifts of atomic transition frequencies 
come from two sources: the finite size of the nu- 
clear charge distribution (the "volume" or "field" shift), 
and the finite mass of the nucleus (see, e.g. J^)- 
The energy shift due to recoil of the nucleus is 
{l/2M)p% = (l/2M)(Spi)2. Furthermore this "mass 
shift" is traditionally divided into the normal mass shift 
(NMS) and the specific mass shift (SMS). The normal 
mass shift is given by the operator (l/2M)Sp^, which 
is easily calculated from the transition frequency. The 
SMS operator is (l/Af)Ei<j(pi • pj) which is difficult to 
evaluate accurately. 

The shift in energy of any transition in an isotope with 
mass number A' with respect to an isotope with mass 
number A can be expressed as 

Siy^'^^ = (fcNMS + fcsMs) ( - +FS{rY''^ , (1) 
where the normal mass shift constant is 

fcNMS = -J^ (2) 
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and (r^) is the mean square nuclear radius. The value 
1823 refers to the ratio of the atomic mass unit to the 
electron mass. 
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In this paper we develop a method for calculating the 
specific mass shift fcgMS for atoms with several valence 
electrons. It is worth noting that in this paper we use 
the convention Su^ = — u^. 

Following our previous work, we are looking for an "all 
orders" method of calculation. Again we have found that 
the finite-field scaling method is very useful in this re- 
spect. The rescaled SMS operator is added to the many- 
particle Hamiltonian 

Hx = Ho + XHsMS ^Ho + xJ2p^- Pr (3) 

The eigenvalue problem for the new Hamiltonian is solved 
for various A, and then we recover the specific mass shift 
constant as 

dE 

Asms = 1™ -tt- (4) 

The operator Q has the same symmetry and structure as 
the initial Hamiltonian Hq (see the Appendix, Ref. [lTI|'). 

We have also calculated the relativistic shift for tran- 
sitions in magnesium using a similar method. These are 
needed in order to measure variation of a Q, H, Q |H IS - 
The dependence of transition frequencies on a can be ex- 
pressed as 

w = + qx, (5) 

where x = (a/ao)^ — 1, ao is the laboratory value of the 
fine-structure constant, and w and wq are the frequen- 
cies of the transition in the absorption clouds and in the 
laboratory, respectively. We vary a directly in codes and 
hence calculate q. This method is described in detail in 
our earlier works (Refs. [Hill El). 

To calculate the energy, we use a combination of con- 
figuration interaction (CI) and many-body perturbation 
theory (MBPT), as was done in Refs. [13, [3. We gen- 
erate a basis set that includes the core and valence or- 
bitals and a number of virtual orbitals. Then we do the 
full configuration interaction calculation in the frozen- 
core approximation. Core-valence correlation effects are 
taken into account using MBPT by altering the Coulomb 
integrals {Q = 1/ \ri — rj\+\pi-pj) in the CI calculation. 
The MBPT operator is separated into two parts: Si and 
S2, which include all one- valence-electron diagrams and 
two-valence-electron diagrams, respectively. The opera- 
tors are calculated to second-order, which was shown to 
be sufficient for single-valence-electron atoms when used 
with a finite-field method pH . 

III. CALCULATION AND RESULTS 

The CI part of the calculation is very similar to the 
large basis set calculation in Ref. We first solved 

the Dirac-Fock equations for the core electrons; we use 
the V^~'^ approximation, so the core includes the or- 
bitals Is 1/2, 2s 1/2, 2pi/2, and 2p^/2- For valence and 



TABLE L Energy of Mg I levels relative to the ground state 
(3s2 ^So). 

Energy (cm~^) 
Level Experiment CI CI -I- Si CI -I- Si + S2 



3s3p 


'6t)0 


21850 


20910 


21676 


21772 


3s3p 


Zpo 


21870 


20930 


21698 


21794 


3s3p 




21911 


20971 


21742 


21837 


3s3p 


Ipo 


35051 


34491 


35474 


35050 


3s4s 


^Sl 


41197 


40406 


41469 


41126 


3s4s 


'So 


43503 


42667 


43744 


43431 


3s3rf 




46403 


45123 


46475 


46306 


3s4p 


3po 


47841 


46919 


48079 


47756 


3s4p 


3po 


47844 


46923 


48082 


47760 


3s4p 




47851 


46929 


48090 


47767 


3s3d 




47957 


46973 


48227 


47880 


3s3ci 




47957 


46973 


48227 


47879 


3s3d 




47957 


46973 


48227 


47879 


3s4p 


Ipo 


49347 


48490 


49672 


49277 


3s4d 




53135 


52041 


53337 


53037 


3s4d 




54192 


53243 


54486 


54111 


3p2 3 


Po 


57813 


56182 


58003 


57706 



virtual orbitals we used a basis formed by diagonaliz- 
ing the Dirac-Fock operator on the basis set of B-splines 
and excluding orbitals with high energy (for a descrip- 
tion of this method as applied in atomic physics, see, 
e.g. Refs. lilSEI). 

The full two-electron CI uses the basis 17 spdf, which 
includes the orbitals 1 — 17si/2; 2 — 17pj, 3 — 17dj, and 
4 — 17/j. It is very close to the saturation of the valence 
CI. The MBPT basis can be larger, since the calculation 
merely involves a summation over the virtual orbitals; we 
have therefore used the basis 32spdfg for this part of the 
calculation, which is essentially complete. 

Table^shows that our ab initio calculation of the spec- 
trum is within 0.4% of the experimental spectrum for all 
considered levels. The relativistic shifts (q-values) are 
presented in Table Table IIIII presents the resulting 
SMS level shift constants, ksus of Eq. (0). In each ta- 
ble we present results of the pure CI calculation (which 
agree with our previous calculation Ref. 0), as well as 
calculations including Si only, and both Si and S2. 

It is worth noting a few points. Firstly, the core- 
valence effects, included using MBPT, make little differ- 
ence to the g-values (less than 10%). This again justifies 
the fact that in previous works for atoms with several va- 
lence electrons these have either been neglected entirely 
or included using a simple fitting procedure based on the 
polarisability of the core 16] . 

The core- valence effects are much more important for 
the SMS calculation. In particular the single-valence- 
electron diagrams (included in Si) can improve accuracy 
drastically in cases where the pure CI method is not very 
good. Although S2 is important for energy calculation, 
it appears to make little difference to fcgMS- This is easily 
understood since the most important two-body diagram 
(the direct diagram, corresponding to the screening of the 
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TABLE II: Relativistic shift of Mg I transitions relative to 
the ground state (3s^ ^So)- 



Upper 


Energy 




q (cm 


Level 


(cm-i) 


CI 


CI + El 


CI + El + Ez 


3s3p ^Po° 


21850 


61 


66 


68 


3s3p ^Pf 


21870 


81 


88 


89 


3s3p ^Pi 


21911 


122 


132 


133 


3s3p ^Pi° 


35051 


86 


94 


94 


3s4s ^Si 


41197 


55 


61 


61 


3s4s % 


43503 


60 


65 


66 


3s3d 


46403 


117 


123 


122 


3s4p 3Po° 


47841 


67 


73 


73 


3s4p ^Pi° 


47844 


70 


76 


77 


3s4p ^Pj 


47851 


77 


83 


84 


3s3d ^Ds 


47957 


79 


85 


86 


3s3d -'Di 


47957 


79 


85 


86 


3s3d ^Di 


47957 


79 


85 


86 


3s4p ipf 


49347 


80 


86 


87 


3s4d ^P>2 


53135 


94 


101 


102 


3s4d ^P>i 


54192 


73 


79 


80 


3p' 'Po 


57813 


198 


214 


214 



TABLE III: Calculations of the specific mass shift constants 
ksMS for Mg I transitions to the ground state (in GHz • amu) . 



Upper 


Energy 




fcsMS (GHz. amu) 


Level 


(cm-^) 


CI 


CI + El 


CI + El + 


3s3p ^Po° 


21850 


-378 


-487 


-492 


3s3p ^Pi° 


21870 


-377 


-486 


-491 


3s3p ^P2° 


21911 


-375 


-485 


-489 


3s3p iPi° 


35051 


231 


120 


134 


3s4s ^Si 


41197 


43 


-59 


-49 


3s4s ^So 


43503 


13 


-94 


-85 


3s3d ^D2 


46403 


-345 


-500 


-477 


3s4p 3Po° 


47841 


-17 


-136 


-126 


3s4p ^Pi° 


47844 


-16 


-136 


-126 


3s4p 3P2° 


47851 


-16 


-136 


-126 


3s3d 


47957 


52 


-87 


-77 


3s3d ^Pi2 


47957 


52 


-87 


-77 


3s3d ^73i 


47957 


52 


-87 


-77 


3s4p iPi° 


49347 


5 


-120 


-108 


3s4d ^Di 


53135 


-100 


-246 


-239 


3s4d ^Di 


54192 


32 


-99 


-88 


3p2 3Po 


57813 


-225 


-469 


-464 



imately 20-30 MHz, which is less than the experimental 
uncertainty in most cases and is of the order of the error 
in our SMS calculations. Furthermore despite numerous 
convergent theoretical calculations of the field shift pa- 
rameter {F of Eq. (^), our knowledge of the field shift 
in magnesium is limited by our lack of knowledge of the 
change in mean-square nuclear radius (5 (r^) ' . 

Also presented in Table IIVI for a theoretical com- 
parison, are the results of Veseth (Ref. 1987) and 
Jonsson et al. (Ref. 1999). Veseth used non- 
relativistic many-body perturbation theory within the al- 
gebraic approximation to calculate the isotope shift to 
third order for some transitions. Jonsson et al. used 
a non-relativistic multiconfiguration Hartree-Fock ap- 
proach, which allowed for both core and valence exci- 
tations in the CI. 

An under-studied transition that is seen in quasar ab- 
sorption spectra is the 2026 A line of Mgl (3s^ 
3s4p ^Pi). From Table Ulll we calculate the isotope shift 
of this hne as 5i/26.24 = 29 50(50) MHz (the error here is 
based on the absence of field shift as well as the incom- 
pleteness of saturation of the basis set used to calculate 

^SMs)- 



IV. CONCLUSION 

We have presented a method for the calculation of the 
isotope-shift in many-electron atoms. It is based on the 
finite-field method, with an energy calculation that com- 
bines CI for the valence electrons and MBPT for the 
core-valence correlations. We have tested the method 
in magnesium, and the agreement was found to be very 
good for all transitions. In particular, for the purposes of 
resolving systematic errors in the search for a-variation, 
and for studies of isotopic evolution of the universe, such 
accuracy is high enough. 

We have also used the method to generate more precise 
values for the relativistic shift (g- values). These were 
found to be within 10% of those found using previous 
methods, as expected. 



electron-electron interaction by the core electrons) makes 
no contribution to the SMS. The exchange diagrams in 
S2 do have an effect, but this is much smaller than the 
one-body contribution. 

In Table IIVI we compare experimental and calcu- 
lated frequency shifts between isotopes ^^Mg and ^''Mg 
(^jy26,24-j^ We compare the SMS part only, which is ex- 
tracted from experiment by subtracting the NMS. We 
have ignored the volume shift for simplicity; it is approx- 
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3s3p ^Po° 
3s3p ^Pi° 
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1173 
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-1409 
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1154 
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